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Abstract 

A study of the effect of spanwise variation in mo- 
mentum on leading edge heat transfer is discussed. 
Numerical and experimental results are presented for 
both a circular leading edge and a 3:1 elliptical lead- 
ing edge. Reynolds numbers in the range of 10,000 
to 240,000 based on leading edge diameter are inves- 
tigated. The surface of the body is held at a constant 
uniform temperature. Numerical and experimental re- 
sults with and without spanwise variations are pre- 
sented. Direct comparison of the two-dimensional re- 
sults, that is, with no spanwise variations, to the ana- 
lytical results of Frossling is very good. The numerical 
calculation, which uses the PARC3D code, solves the 
three-dimensional Navier-Stokes equations, assuming 
steady laminar flow on the leading edge region. Ex- 
perimentally, increases in the spanwise- averaged heat 
transfer coefficient as high as 50 percent above the 
two-dimensional value were observed. Numerically, 
the heat transfer coefficient was seen to increase by as 
much as 25 percent. In general, under the same flow 
conditions, the circular leading edge produced a higher 
heat transfer rate than the elliptical leading edge. As 
a percentage of the respective two-dimensional values, 
the circular and elliptical leading edges showed sim- 
ilar sensitivity to spanwise variations in momentum. 
By equating the root mean square of the amplitude 
of the spanwise variation in momentum to the turbu- 
lence intensity, a qualitative comparison between the 
present work and turbulent results was possible. It 
is shown that increases in leading edge heat transfer 
due to spanwise variations in freestream momentum 
are comparable to those due to freestream turbulence. 

Nomenclature 

(pV) momentum 

a area 

A peak to peak momentum difference 

normalized by the characteristic momentum 
6 blockage, defined as ratio of maximum 

* Member AIAA 


body thickness to tunnel height 
D leading edge diameter (twice the leading 

edge radius of curvature) 

Fr Frossling number ( Fr = Nu/y/Rto) 

h heat transfer coefficient ( h = q/{T w — T au) 

k thermal conductivity of fluid 

M Mach number 

Nu Nusselt number (Nu = hD/k ) 

p pressure 

P refers to the position 2 diameters 

directly above the stagnation point 
Pr Prandtl number 

q surface heat flux 

Re d Reynolds number based on leading 

edge diameter (Rep = pVD/p) 
s distance along body surface 

T temperature 

Tu turbulence intensity (also Tu = ApjyJ 2) 

V speed 

xo distance between wire array and 

body leading edge 

Greek Symbols 

7 ratio of specific heats 

p kinematic viscosity of fluid 

A wavelength of spanwise variation 

p density of fluid 

Subscripts and Superscripts 

()o two-dimensional results (i.e., A = 0) 

Oau; adiabatic wall 

() c /j characteristic quantity 

Qedge conditions at edge of boundary layer 

Oexit conditions at the exit 

() tn conditions at the inlet 

() 5t conditions at stagnation point 

()t total conditions 

() w conditions on body (wall) 

() spanwise-averaged quantity 


This paper is declared a work of the U.S. Government and 
is not subject to copyright protection in the United States 
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I Introduction 

Accurate prediction of the heat transfer rates on 
gas turbine blades is very important to the design of ef- 
ficient and reliable engines. Typically, the leading edge 
region endures heat transfer rates comparable to the 
peak value for the entire blade. The leading edge heat, 
and analogously mass transfer, is known to be sensi- 
tive to even low levels of freestream turbulence, and to 
small amplitude spanwise variations in the freestream 
momentum [1, 2, 3, 4, 5, 6]. However, the mechanisms 
responsible for this sensitivity are not well understood. 
Attempts to explain this sensitivity in terms of a pul- 
sating [7] or oscillating [8] freestream have not pro- 
duced increases comparable to those seen experimen- 
tally. Rigby and VanFossen [9] studied, numerically, 
the effect of spanwise variation of freestream momen- 
tum on the heat transfer to a cylindrical leading edge. 
They found that stretching of the freestream vorticity 
introduced by the spanwise variation created highly 
three-dimensional flow in the leading edge region. This 
vortex structure near the leading edge was determined 
to cause large increases in leading edge heat transfer. 

It is hoped that a better understanding of the 
mechanisms which allow low levels of freestream turbu- 
lence to cause large increases in stagnation heat trans- 
fer can result from the study of a particular isolated 
phenomenon. The present study addresses the effect 
of leading edge contour on heat transfer to a body in 
crossflow. Results with and without spanwise varia- 
tions in the freestream momentum are discussed. Nu- 
merical and experimental results for a cylindrical lead- 
ing edge and a 3:1 elliptical leading edge are presented. 
In addition, a specific application of the analysis of 
Frossling [10] is used to predict the two-dimensional 
heat transfer to each of the bodies investigated. 


II Statement of Problem 

Consider the problem of a body with a blunt lead- 
ing edge. Figure 1 shows, schematically, the geometry 
of the problem (not drawn to scale). The body is uni- 
form in the spanwise direction ana extends far down- 
stream of the leading edge. Let D define the leading 
edge diameter at the geometric stagnation point. This 
diameter is defined to be twice the leading edge radius 
of curvature. The flow is confined on four sides by 
walls which define the test section. Note that the test 
section does not vary in width or height in the stream- 
wise direction. The surface of the body is held at a 
uniform temperature. In the freestream a spanwise 
variation in momentum is introduced, characterized 
by a wavelength, A, and an amplitude, A. The line 
from P to PP in figure 1 is at the same streamwise lo- 
cation as the leading edge, and is two diameters above 
the centerline of the tunnel. The amplitude of the 
spanwise variation of momentum on this line P — P P 
will be referred to as Ap, Experimentally, this span- 
wise variation is produced by a vertical array of wares 
located upstream of the body. Numerically, the span- 
wise variation in momentum is specified as a sinusoidal 
variation at the upstream boundary. 

The two geometries to be considered are a circular 
leading edge and a 3:1 elliptical leading edge. Figure 1 
emphasizes the circular leading; the 3:1 elliptical lead- 



Figure 1: General problem (not drawn to scale). 

ing edge is shown as a dashed line for reference. In 
each case, the leading edge radius of curvature is the 
same. Thus, for the circular leading edge, the plate 
thickness is D. For the 3:1 elliptical leading edge, the 
plate thickness is 3D. 

Ill Numerical Method 

To solve this problem numerically, the PARC3D 
code[ll] was used. PARC3D uses the Beam- Warming 
approximate factorization algorithm and a time- 
marching implicit scheme which solves the Navier- 
Stokes equations in strong conservation form. Some 
of the basis assumptions used in the PARC3D code 
are that the fluid behaves as a perfect gas, that the 
Prandtl number is constant, and that the temperature 
dependence of the kinematic viscosity is approximated 
by Sutherland’s law [14, p. 328]. PARC3D uses second 
order accurate central-differencing for spatial deriva- 
tives. For stability, fourth difference artificial dissipa- 
tion is added. The coefficient of the artificial dissi- 
pation varies both spatially and directionally, which 
produces a stable and accurate solution. This form of 
the artificial dissipation is especially important when 
computing flov^s with vorticity in the freestream, since 
excessive amounts of artificial dissipation can acceler- 
ate the decay of vorticity. In the present study, care 
has been taken to ensure that decay of freestream vor- 
ticity is dominated by the physical phenomena, and 
not by the numerical method. This is done by com- 
paring the rate of decay of vorticity upstream of the 
body to an analytical result which neglects the pres- 
ence of the body [9]. 

The boundary conditions used are as follows. On 
the body surface, a no-slip condition wuth uniform 
temperature is applied. At the' inlet, a spanwise to 
tal pressure variation is imposed so as to produce a 
sinusoidal momentum profile, while inlet total tem- 
perature is uniform. At the exit, static pressure is 
specified, while other flow quantities are extrapolated. 
Symmetry is imposed in the spanwise direction at half 
wavelength intervals. Thus, only one half wavelength 
in the spanwise direction is computed. Symmetry is 
also imposed above and below the body by assuming 
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it is at zero angle of attack. Thus, only the top half of 
the flow field is computed. To take into account the 
upper wall of the test section, an inviscid wall bound- 
ary condition was used. That is, the flow was required 
to remain tangent to the upper surface with no flow 
passing through it. No attempt was made to resolve 
the boundary layer on this upper surface. 

It has been shown in a previous study [9] that the 
PARC3D code can produce accurate solutions to this 
type of problem. Grids used in the present work are 
even finer than those that had been shown to be ade- 
quate in Rigby and VanFossen [9]. For each case, the 
residual was required to drop three orders of magni- 
tude from the initial level. The residual is defined to be 
the sum of the L2 norms for each of the five equations, 
divided by the number of grid points. In addition, the 
value of the Frossling number at the leading edge was 
monitored. No solution was deemed to have converged 
until the per iteration change of this Frossling number 
fell below 10“ 6 . 

For the present application, PARC3D was set up 
to solve the steady laminar three-dimensional flow 
equations in the leading edge region, the area of in- 
terest in this problem. In general, a turbulent solution 
is required to keep stable any separation which might 
occur near the junction between leading edge and flat 
plate. Initiating a turbulent calculation near the junc- 
tion, while greatly improving the convergence rate, has 
very little effect on the flow solution in the upstream 
laminar region. Thus downstream of the leading edge, 
the algebraic Baldwin-Lomax [12] turbulence model is 
used. 


IV Experimental Method 

Wind Tunnel 

All tests were conducted in a wind tunnel which 
was connected to the laboratory exhaust system. This 
wind tunnel is described in more detail in reference [2], 
While operating, room air flowed into a 73.7 cm wide 
by 68.6 cm high rectangular bell mouth past tur- 
bulence damping screens and soda straws and then 
through a 4.85:1 contraction. The dimensions of the 
test section were 15.2 cm wide by 68.6 cm high. The 
maximum velocity obtainable in the test section was 
46 m/s and the clear tunnel turbulence intensity was 
less than 0.5 percent at all flow rates. Mass flow rates 
were measured with an ASME standard sharp edge 
orifice in the exhaust line; a calibration between mass 
flow and test section velocity measured in the clear 
tunnel with a Pitot tube was used to determine test 
section velocities during runs. For mass flow rates be- 
low the capability of the orifice, a hot wire was used to 
measure test section velocities. Inlet total temperature 
was measured with four thermocouples just outside the 
inlet bell mouth. 

To establish spanwise variations in momentum, a 
vertical array of wires perpendicular to the span was 
placed upstream of the leading edge. The wire array 
was made up of 0.51 mm diameter wires spaced 6.3 
mm apart. The amplitude of the spanwise variation 
reaching the leading edge was varied by changing the 
streamwise location of the array. The array was posi- 
tioned from 15.6 cm to 30.5 cm upstream of the body 


leading edge. The position of the wire array relative 
to the leading edge of the body will be referred to as 
xq . Measurements were also taken with the wire ar- 
ray removed, producing what will be referred to as the 
two-dimensional results. 


Heat Transfer Models 

Spanwise- averaged heat transfer coefficients were 
measured on a circular leading edge with a 3.30 cm ra- 
dius and on a 3:1 elliptical leading edge. The elliptical 
leading edge design had a semi-major axis of 29.72 cm 
and a semi-minor axis of 9.91 cm. Thus the elliptical 
leading edge had a radius of curvature at the stagna- 
tion point of 3.30 cm, which matches the radius of the 
circular leading edge. Both models spanned the 15.2 
cm wide test section. 

Heat transfer coefficients were measured on the 
surface of the models with electrically heated alu- 
minum strips. Each strip was 6.60 cm long in the 
spanwise direction by 4.76 mm wide in the flow direc- 
tion and 3.18 mm deep. A thin film electric heater 
encapsulated between two layers of polymide plas- 
tic was fastened to the back of each aluminum strip 
with a pressure sensitive adhesive; a thermocouple 
was then embedded in a groove on the side of each 
strip. The aluminum heat flux gauges were embedded 
in the surface of each model with a center- to- center 
spacing along the surface of about 5.0 mm. The 0.25 
mm wide gap between each heat flux gauge was filled 
with epoxy. A guard heater was also used behind the 
surface-mounted heat flux gauges to minimize heat lost 
to the interior of the model. With the desired flow 
established, electric current was applied to each gauge 
individually to keep them all at the same temperature. 
Heat transfer coefficients for each strip were then cal- 
culated from the equation 


h = 


QEI Qrad ~ Qend Qgap 
dgagt (Tgage ~~ T a w ) 


( 1 ) 


where qei is the electric power dissipated in the foil 
heater, q ra d is the heat lost by radiation from the gauge 
surface, q en d is the heat lost to the epoxy from the un- 
guarded ends of the aluminum strip, and q gap is the 
heat lost to the epoxy gap between gauges. a gage is 
the area of the gauge surface exposed to the air, T gage 
is the assumed uniform temperature of the gage mea- 
sured by the thermocouple, and T aw is the adiabatic 
wall temperature. The radiation heat loss was esti- 
mated using an emissivity of 0.05 for the aluminum. 
The conduction losses q en d and q gap were estimated 
from heat flow in a rectangle (the epoxy) with one side 
at a constant temperature (the edge of the aluminum), 
an adjacent side with convection to air, and the other 
two sides assumed to be insulated [15, p. 167]. 


V Results 

Two-Dimensional Results 

General Remarks 

In this section, cases with no spanwise vari- 
ation are discussed. One would expect from the 
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Frossling analysis [10] that, for two-dimensional flow, 
the Frossling number results should be nearly indepen- 
dent of Reynolds number. The Frossling number is a 
dimensionless measure of the heat transfer defined by 


Nu 

\/ReD 



J iP v )^D 

V Vch 


( 2 ) 


These two-dimensional results serve two useful pur- 
poses. First, they demonstrate the effect of body shape 
on the heat transfer. Second, the close agreement be- 
tween experiment, numerical calculation, and analysis 
lends confidence to the three-dimensional results. 

When defining dimensionless quantities such as 
Reynolds and Frossling numbers, one must decide on 
the characteristic values of various parameters. Ef- 
fects such as variable properties, blockage, and com- 
pressibility must be addressed. Thus values for these 
parameters were chosen to maximize the applicability 
of the results by minimizing the variation in Frossling 
number from case to case. 

Both the experimental and numerical results expe- 
rience the effects of blockage due to surrounding walls. 
For internal flow, the choice for characteristic momen- 
tum is not immediately apparent. One might choose 
the inlet value, the exit value, or some other value. For 
the present study, a simple average between the inlet 
value and the value in the most blocked region seemed 
to best account for the effects of blockade. Since the 
test section geometry does not vary with streamwise 
distance, the most blocked streamwise position occurs 
at the thickest part of the body. Thus the blockage 
can be defined as follows 

£ _ max, body thickness ^ 

tunnel height 


In terms of the inlet momentum, (pV) C h can be written 
as 

(pV)ch = ( pV)in ^ W 

As mentioned above, this is a simple average between 
the inlet momentum and the momentum which would 
occur at the most blocked region due to conservation 
of mass. 

The effect of variable properties appears to be 
best accounted for by evaluating the viscosity, p, and 
the thermal conductivity, k, at the freestream total 
temperature. Evaluating these properties at a film 
temperature, defined as a simple average between the 
wall and the total temperatures, was also investi- 
gated. It was found that using the film temperature, 
as opposed to using the total temperature, caused the 
Frossling number to have a much stronger dependence 
on temperature difference. This observation was most 
strongly demonstrated by the numerical results, where 
cold and hot wall results could be compared. Wall tem- 
peratures in the range 0.792 T t < T w < 1.1887* were 
calculated. 

It is informative to note the progression of ideas 
which took place. Preliminary hot wall calculations 
without blockage were done for the 3:1 elliptical lead- 
ing edge. Using the inlet momentum for the Reynolds 


number and the freestream total temperature to eval- 
uate properties, it was observed that experimental val- 
ues of Frossling number were higher than the numer- 
ical ones, which agreed well with the Frossling anal- 
ysis. In an attempt to produce better agreement, 
the properties were then evaluated at the film tem- 
perature, which caused the numerical and experimen- 
tal Frossling numbers to be reduced producing better 
agreement between experiment and analysis, but worse 
agreement between numerical and analytical results. 
Note that in the analytical results, constant proper- 
ties are assumed. Experimentally, wall temperatures 
below the freestream value could not be produced. 
However, several numerical several runs with nega- 
tive temperature difference were calculated. For these 
cases, using film temperature instead of freestream to- 
tal temperature resulted in Frossling numbers higher 
than the analytical result. In fact, relative to using the 
total temperature, the correction from using the film 
temperature was in the wrong direction. 

Finally, it was realized that using the total temper- 
ature to evaluate viscosity and thermal conductivity 
produced nearly the same result for both positive and 
negative temperature differences. The values were, 
however, consistently higher than those measured ex- 
perimentally. At this point, the correction due to 
blockage was introduced, yielding the desired minimal 
variation in Frossling number and good agreement for 
experimental, numerical, and analytical results. Recall 
that this statement refers to a comparison between nu- 
merical results without blockage and the experiment 
which had blockage of approximately 0.3 for the 3:1 
elliptical leading edge. At this point it was decided 
that all numerical calculations would be run with the 
effects of blockage included. With the blockage nearly 
matched between the numerical and experimental re- 
sults, agreement is good with or without the correction 
for blockage. However, by leaving the effect of blockage 
in the definition for the characteristic momentum, it 
is believed that the results presented will more closely 
correspond to results for unblocked flow. 

Eckert and Drake [13, pp. 263-266] suggest that 
the most appropriate characteristic temperature dif- 
ference for the evaluation of the heat transfer coeffi- 
cient h is the difference between the wall temperature 
and the adiabatic wall temperature. They suggest the 
adiabatic wall temperature can be approximated by 

T aw = V/MT, - T edge ) + T edge (5) 

The edge temperature is deduced from the surface 
pressure distribution obtained from the numerical re- 
sults. Using the isentropic relations, which are valid 
outside the boundary layer, and assuming that the 
pressure does not vary across the boundary layer, the 
edge temperature can be written as 

T edge = T t ( E ^ L ) (6) 

Numerical 

Calculations were done using the PARC3D 
Navier-Stokes code. Boundary conditions and the grid 
geometry are specified to mimic the experimental fa- 
cilities. Figure 2 shows the grid for the 3:1 elliptical 
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leading edge. As can be seen in figure 2, the grid is one 
half of a C-type grid. Symmetry was imposed across 
the stagnation point by assuming a zero angle of at- 
tack. The two-dimensional calculations were done on 
a 48x61x5 grid. Five grid points are required in the 
spanwise direction because of the three-dimensional 
nature of the PARC3D code. These two-dimensional 
calculations were done with a Reynolds number rang- 
ing from 10,000 to 200,000, and with a wall tempera- 
ture from 0.792T* to 1.188T t . Additionally, the inlet 
Mach number was varied from 0.1 to 0.2. 



Figure 2: 48x61 Grid. 

Recall that the numerical calculation assumes a 
constant value for Prandtl number throughout the 
flowfield. The default value of the Prandtl number in 
the PARC3D code is 0.72, and this is the value which 
was used for all the calculations. For the present ex- 
perimental results, a Prandtl number of 0.71 is more 
appropriate. Rather than rerun calculations which had 
already been done, an adjustment to the Prandtl num- 
ber of 0.71 is used [16]. All numerical values of the 
Frossling number are adjusted by assuming that 

^ * Fr ° " (sir <7) 

Thus, all results presented in this discussion are ad- 
justed to a Prandtl number of 0.71. 

Table 1 provides a summary of the cases run with 
PARC3D and the resulting Frossling number at the 
stagnation point. The minimal variation in the stagna- 
tion point value of the Frossling number indicates that 
appropriate decisions have been made for the char- 
acteristic quantities. Figure 3 shows the variation of 
the Frossling number with streamwise distance along 
the body, for the circular and the 3:1 elliptical leading 
edge. Note that all of the cases listed in table 1 are 
plotted in figure 3, which also shows the results of a 
specific application of the Frossling analysis [10] de- 
scribed in the appendix. The reader is strongly urged 
to refer to the appendix to see how the Frossling anal- 
ysis can be applied for specific body contours. Agree- 
ment between these numerical results and the Frossling 
analysis is good. 

The minimal variation of Frossling number with 
increasing Reynolds number indicates that grid reso- 
lution normal to the body surface is sufficient for these 
cases. For the 3:1 elliptical leading edge, there were 17 
grid points in the boundary layer at Rto = 10, 100, 



Ren 

'Ll. 

Tt 

Mch 

Fr st 

elliptical 

10,100 

1.075 

0.18 

0.781 

elliptical 

20,300 

1.075 

0.18 

0.766 

elliptical 

51,000 

1.075 

0.18 

0.775 

elliptical 

101,900 

1.075 

0.18 

0.777 

elliptical 

203,800 

1.075 

0.18 

0.778 

circle 

23,500 

1.074 

0.21 

0.939 

circle 

23,500 

1.089 

0.21 

0.939 

circle 

23,500 

0.891 

0.21 

0.935 

circle 

23,500 

0.911 

0.21 

0.936 

circle 

23,500 

0.792 

0.21 

0.932 

circle 

23,500 

1.188 

0.21 

0.942 

circle 

24,600 

1.082 

0.11 

0.939 

circle 

24,600 

0.918 

0.11 

0.934 


Table 1: Summary of two-dimensional calculations. 



Figure 3: Frossling number vs distance, two- 

dimensional calculations listed in table 1. 

and 11 grid points in the boundary layer for Rep — 

203.700. As will be shown below, agreement with ex- 
periment is good, which indicates that the streamwise 
grid resolution is also sufficient for these body shapes. 

Experimental 

Experimentally, clear tunnel results were ob- 
tained for Reynolds numbers ranging from 77,500 to 

232.700. Wall temperatures on the average of 1.076 
times the freestream total temperature are presented. 
Because the freestream total temperature was the am- 
bient room temperature, approximately 18° Celsius, 
the temperature difference was approximately 22° Cel- 
sius above the freestream total. Table 2 shows a sum- 
mary of the clear tunnel results. The stagnation value 
of Frossling number is seen to vary little with changes 
in Reynolds number and temperature difference. Fig- 
ure 4 shows the variation of Frossling number with 
distance along the body. This figure shows the re- 
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suits from each of the cases shown in table 2. Also 
on this figure are the results of the PARC3D calcula- 
tion. Agreement is very good over the measured re- 
gion. These results are also in good agreement with 
similar experimental results for a circular cylinder [14, 
p. 305], and for a 3:1 ellipse [17] in crossflow. 



Rto 

T m 

T t 

Fr, t 

elliptical 

77,500 

1.074 

0.750 

elliptical 

108,000 

1.074 

0.750 

elliptical 

150,900 

1.074 

0.770 

elliptical 

188,400 

1.073 

0.777 

elliptical 

232,700 

1.074 

0.759 

circle 

189,600 

1.078 

0.923 

circle 

117,400 

1.077 

0.927 

circle 

151,100 

1.078 

0.934 

circle 

189,600 

1.078 

0.923 


Table 2: Summary of two-dimensional experimental 
results. 



Figure 4: Frossling number vs distance, two- 

dimensional measurements listed in table 2. 

Looking at both figures 3 and 4 it can be seen 
that agreement between the numerical, experimental, 
and analytical results is very good. This agreement 
reinforces the confidence in each of the methods. Also 
notice that the Frossling number for the 3:1 elliptical 
leading edge is roughly 20% lower than for the circular 
leading edge. The lower value is a direct result of the 
lower gradient in edge momentum at the stagnation 
point for the 3:1 elliptical leading edge. The reader 
is strongly urged to see the appendix for a complete 
explanation of this statement. 


Three-Dimensional Results 

Numerical 

Numerical results were produced for Reynolds 
numbers from 22,000 to 100,000. A sinusoidal span- 
wise profile in momentum was imposed 4.5 D upstream 
of the leading edge. The wavelength of the spanwise 
profile was 0.1Z) for all cases, which matches the wire 
spacing in the experiment. Cases were run with am- 
plitudes in the range of 1 to 15% of the freestream mo- 
mentum. The numerical grid used for the three dimen- 
sional runs matched that used for the two-dimensional 
calculations, except 32 grid points were used in the 
spanwise direction. The x - y grid was simply stacked 
in the z direction, equally spaced over one-half wave- 
length. As in the two-dimensional calculations, sym- 
metry was imposed across the stagnation point. In 
addition, symmetry was imposed at every half wave- 
length in the spanwise direction. The Mach number at 
the inlet for all cases was approximately 0.15, which 
was about as low as the PARC3D code could go and 
still produce reasonable convergence times. No at- 
tempt was made to match the Mach number in the 
experiment since, at these low values, it is not ex- 
pected to affect the results. Figures 5 and 6 show 
the variation of spanwise- averaged Frossling number 
along the body surface, for the circular and ellipti- 
cal leading edges, respectively. In each figure, results 
for several amplitudes are shown as well as the result 
from the Frossling analysis. Recall that the amplitude 
Ap noted in these figures refers to the amplitude of 
the momentum variation at position P shown in fig- 
ure 1. The Frossling number was increased above the 
two-dimensional value for every case, with increases 
as high as 25% observed. As a percentage of the 
corresponding two-dimensional value, the circular and 
the elliptical leading edges show similar increases in 
Frossling number for equivalent freestream conditions. 
Notice that heat transfer is increased quite uniformly 
for s/D < 0.5. That is, each curve appears similar 
to the two-dimensional results, but is shifted up by a 
certain amount. Note that the spike in figure 5 occurs 
near the junction between the circle and the flat plate 
portion of the body. This is also the position where 
the turbulent solution is initiated. However, no spike 
is seen in the 3:1 elliptical leading edge results since the 
junction between the ellipse and the flat plate occurs 
at $/D & 4. 

Experimental 

Measurements were taken for Reynolds number in 
the range of 13,000 to 240,000. The magnitude of the 
spanwise variation in momentum reaching the body 
was varied by positioning the array of wires at var- 
ious streamwise locations. Figures 7 and 8 show the 
variation of spanwise- averaged Frossling number along 
the body surface for the circular leading edge and the 
3:1 elliptical leading edge, respectively. Increases in 
the spanwise-averaged Frossling number by as much as 
50% above the two-dimensional value were observed. 
Introduction of the wires is seen to cause a uniform in- 
crease in the heat transfer over most of the measured 
region. These increases are substantial, and demon- 
strate the acute sensitivity that the leading edge heat 
transfer has to small spanwise variations in freestream 
momentum . 





Figure 5: Frossling number vs distance, 3D calcula- 
tions, circular. 



Figure 6: Frossling number vs distance, 3D calcula- 
tions, elliptical. 


With the 3:1 elliptical leading edge model in place 
and the array of wires located 2.79 D upstream of the 
leading edge, spanwise hot-wire traverses were used 
to measure the momentum profiles. Figure 9 shows 
the measured amplitude at position P, expressed as a 
percentage of the reference momentum. Position P is 
located at the same streamwise location as the leading 
edge, and approximately two leading edge diameters 
above the tunnel centerline, see figure 1. Notice that 
the amplitude is relatively high at low Reynolds num- 
ber, then drops off and settles on 4% over much of 
the Reynolds number range, perhaps due to mixing 
of wire wakes. Figure 10 shows the hot-wire trace for 
one of the largest amplitudes. In this trace, one can 
see a dominant spanwise variation with spacing equal 
to the spacing in the array of wires. In addition, there 
is an unsteady component resulting from low levels 
of freestream turbulence and unsteadiness in the wire 
wakes. The amplitude for each trace was arrived at by 
simply looking at the trace and estimating the average 
peak-to-peak variation. 



Figure 7: Frossling number vs distance, 3D measured, 
circular, xq = 2.36D. 



Figure 8: Frossling number vs distance, 3D measured, 
elliptical, xo = 2.79D. 

Discussion 

In this section a more direct comparison between 
measured and calculated heat transfer is presented. It 
has been demonstrated above that both the numer- 
ical and experimental results show a rather uniform 
increase in Frossling number over the leading edge re- 
gion. Because of this uniform increase, attention will 
be focussed on the value at the stagnation point. In 
addition, each Frossling number is normalized by the 
appropriate value from the Frossling analysis described 
in the appendix. For the circular leading edge, the 
value is 0.9411. For the elliptical leading edge, the 
value is 0.7691. 

It is informative to make a qualitative comparison 
between turbulent results and the present results. To 
do this, the root mean square value of the amplitude 
in the present study is equated with the turbulence 
intensity. Assuming a sinusoidal spanwise variation 
in momentum, the root mean square value is simply 
Ap/y/2. As a comparison to turbulent data in general, 
the correlation of Kestin and Wood [19] is used. While 
there is much scatter in the heat transfer data in the 
literature, this correlation does give a qualitative indi- 
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Figure 9: Experimental Ap versus Reynolds number, 
elliptical, xo = 2.791?. 



Figure 10: Sample hot-wire trace, elliptical, xq = 
2.79D, Re D = 13,000. 


cation of the trend observed in turbulent flow over a 
circular cylinder. The correlation, normalized by the 
stagnation point value, is expressed as 

3.48 fTuy/ReE\ 3.99 f Tu^/ReE ^Y 
Fr ~ l + 0.945 V 100 ) 0945 V 10° / 

( 8 ) 

Figure 11 shows some of the heat transfer results 
plotted as a function of the quantity Tuy/Rep. Also 
shown in this figure is the correlation of Rest in and 
Wood, equation 8. The experimental results shown in 
figure 11 are for the 3:1 elliptical leading edge, with 
the array of wires located 2.361? upstream. This is 
the only set of measurements for which amplitude lev- 
els are available. Notice that the present experimen- 
tal and the numerical results have a trend similar to 
the correlation. The exception is that the numerical 
results seem to level off after an increase of roughly 
25% is achieved. The numerical results indi cate a very 
strong sensitivity at low values of Tu y/Rep . However, 
it would appear that for Tuy/Rep > 5, additional 
flow physics, such as unsteadiness, may be required 
for direct comparison to experiment. In spite of this, 
a substantial part of the increase in heat transfer does 
seem to be accounted for in t he nu merical results, even 
at the higher values of Tuy/Rep . Another important 
observation about the numerical results shown in fig- 
ure 11 is that the circular and elliptical results plot 


nearly on top of one another when normalized by their 
respective two-dimensional values. As will be shown 
below, the experimental results also exhibit this be- 
havior. If this result is true in general, then know- 
ing the reaction of one body shape to a certain flow 
could provide information about other body shapes, 
which could expand the applicability of many existing 
results. Experimentally, the smallest increase in the 
Frossling number corresponds to the lowest Reynolds 
number, while the largest increase occurrs at the high- 
est Reynolds number. 



Figure 11: Stagnation Frossling number vs Tuy/Rep ■ 

Figure 12 shows the same data as in figure 11, 
except plotted as a function of Ap } the momentum 
amplitude at position P. As stated above, position 
P is located at the same streamwise location as the 
leading edge of the body, and approximately two lead- 
ing edge diameters above the tunnel centerline. This 
plot demonstrates the very sharp rise in the numer- 
ical Frossling number for amplitudes as low as 1 to 
2% of the reference momentum. The numerical re- 
sults exhibit a much stronger sensitivity to Ap than 
to Rep. Due to the number of points on figure 12, 
Reynolds numbers are not indicated for each point. 
However, the relative insensitivity to Rep can be seen 
by focussing one’s attention on the six points clustered 
above Ap = 2%. These six points range in Reynolds 
number from 22,000 to 100,000. In fact, the largest 
increase at Ap = 2% occurred for the lowest Reynolds 
number. On the other hand, the experimental results 
exhibit quite a strong effect as Reynolds is varied. Re- 
ferring once again to figure 12, one sees that as the 
Reynolds number increases, the amplitude in the ex- 
perimental results is fairly constant once it drops off, 
and yet the heat transfer keeps rising. It has previ- 
ously been shown [9] that flowfields such as the one 
studied here may nave a spanwise periodic array of 
vortex structures situated just in front of the stagna- 
tion region. These vortex structures are observed in 
the present calculations. It is observed that as the vor- 
tex structure becomes stronger, its position is farther 
upstream, thus diminishing its ability to increase the 
heat transfer. While it is observed numerically that 
the vortex structure becomes stronger with increasing 
Reynolds number, the previous statement may explain 
a relative insensitivity of the heat transfer to changes 
in Reynolds number for the numerical results. 
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Ap, % 

Figure 12: Stagnation Frossling number vs Ap. 

Figures 13 and 14 show the span wise- averaged 
stagnation point Frossling number as a function of the 
square root of the Reynolds number. Also shown in 
these figures is the correlation of Kestin and Wood, 
equation 8. As mentioned above, for the sake of com- 
parison, Tu is equated with Ap/y/ 2. Since the ampli- 
tude for the experiment is around 4% over most of the 
Reynolds number range, as shown in figure 9, a value 
of Tu is chosen to be Tu = (0.04)/\/2 = 0.028. Each 
result in figures 13 and 14 is normalized by the respec- 
tive two-dimensional value. The experimental results, 
figure 14, have a kink near Re& = 40,000 (i.e. 200 2 ), 
then continue to increase as the Reynolds number in- 
creases. As was seen in figure 9, the momentum varia- 
tion reaching the body initially drops with increasing 
Reynolds number, then levels off. This drop in ampli- 
tude at the lower Reynolds numbers could explain the 
leveling off seen in the experimental Frossling number 
at the lower Reynolds numbers. It must also be noted 
that at the lowest flow rate, the wire wakes appeared 
to be steady. That is, there was no vortex shedding 
in the wire wakes. As Reynolds number was increased 
above 20,000 the wakes became unsteady. Thus, it can 
not be entirely ruled out that wire wake unsteadiness 
may be having some effect on these results. 

In general, for the numerical results, a larger am- 
plitude produced a larger increase in heat transfer for 
a given Reynolds number. However, increases of more 
than about 25% of the corresponding two-dimensional 
value does not appear possible for the present calcu- 
lations. Perhaps additional flow physics, such as un- 
steadiness, must be included at higher Reynolds num- 
bers. In the calculation of horseshoe vortex flow, Vis- 
bal [18] found that a time-periodic solution was pos- 
sible, even with steady boundary conditions. The 
present flow situation is similar to the horseshoe vortex 
problem in some ways. However, the PARC3D code 
was not run time-accurate, so all solutions presented 
represent steady flow. Perhaps for an high enough 
Reynolds number with an high enough amplitude, an 
unsteady flow for the present geometry could result. 
The levels required for unsteadiness to occur could not 
be predicted at this time. 

One final and very interesting result shown in fig- 
ures 13 and 14 is that, when normalized by their re- 
spective two-dimensional values, the circular leading 



Figure 13: Stagnation Frossling number vs y/Rep, 
PARC3D. 



Figure 14: Stagnation Frossling number vs \/ Re d. 
Measured. 


edge and elliptical leading edge show similar increases 
in heat transfer for similar flow conditions. This simi- 
lar increase in heat transfer is seen in both numerical 
and experimental results. If this result is generally 
true, then knowing the reaction of one body shape to 
a certain flow situation could be used to predict the 
reaction of a second body to the same flow situation. 


VI Conclusions 

It has been shown numerically and experimentally 
that spanwise variation in the freestream momentum 
is a very strong mechanism for increasing leading edge 
heat transfer. For small spanwise variations or for low 
Reynolds numbers, the numerical results showed in- 
creases in heat transfer much like those seen exper- 
imentally. However, these three-dimensional steady 
laminar numerical results were only able to account 
for at most a 25% increase in leading edge heat trans- 
fer. Experimental results showed increases as high as 
50% above the two-dimensional value. It is felt that 
increases above 25% may be due to unsteady effects, 
such as unsteady wire wakes or time- varying vortex 
structures near the leading edge. Equating the root 
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mean square of the amplitude in the present study 
to the turbulence intensity of other results in the lit- 
erature showed qualitative agreement. That is, the 
leading edge heat transfer showed a similar sensitiv- 
ity to small steady spanwise variations in freestream 
momentum, as it did to freestream turbulence. 

It was also seen that two bodies with equivalent 
flow conditions in the freestream experience similar in- 
creases in heat transfer, relative to their respective 
two-dimensional values. These similar increases in 
heat transfer indicate that heat transfer results for 
a given body, when normalized by the laminar two- 
dimensional value, could be used to predict the heat 
transfer to a different body shape under similar flow 
conditions. 

Two-dimensional results were obtained numeri- 
cally, experimentally, and analytically. Agreement was 
good for all these results, providing confidence in each 
of the methods. For equivalent flow conditions, the 
heat transfer rate to the 3:1 elliptical leading edge was 
for all cases approximately 20% lower than that of the 
circular leading edge. This 20% difference is seen to be 
a consequence of the lower gradient in edge velocity at 
the leading edge for the 3:1 elliptical leading edge. A 
specific application of the analysis of Frossling is seen 
to account completely for the differences due to body 
shape. 

One might consider freestream turbulence to be 
made up of unsteadiness and instantaneous spatial 
variations in the freestream. By studying the effect of 
one particular isolated phenomenon, it is hoped that 
a better understanding of the mechanisms which al- 
low low levels of turbulence to cause large increases in 
stagnation point heat transfer will result. The work of 
other researchers, cited herein, suggests that the effect 
of oscillating or pulsating flow on the stagnation heat 
transfer is quite small. By contrast, the present study 
shows that another phenomenon which would exist at 
any instant in every turbulent flow, namely spanwise 
variations, can cause large increases in leading edge 
heat transfer. Thus further study of this mechanism 
appears warranted, in particular, the combination of 
spanwise variations and unsteadiness. 
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Appendix 

Frossling Analysis 

The analysis of Frossling [ 10 ] is a boundary layer 
approximation for laminar stagnation point flow, as- 
suming incompressibility and constant properties. The 
reader is referred to the literature for the details of the 
analysis; the results are as follows. Let the velocity 
at the edge of the boundary layer, normalized by the 
freestream value, be approximated by the following 

W = + %(£)’ + %(£)' (A.l) 


For application of the Frossling analysis presented 
above, a representative case for each body shape was 
chosen. The edge momentum for each body was de- 
duced from the surface pressure in the numerical cal- 
culation. Using the isentropic relations, and assuming 
constant pressure across the boundary layer, the edge 
momentum can be written as 

{pV ) tigt = Jiptxit h (P‘ - P*«») ( A - 8 ) 

A least squares regression to a fifth order polynomial 
is then done in the range 0 < $/ D < 1. The result for 
the circular leading edge is 


where s is the distance measured along the surface, D 
is the leading edge diameter, and Vi, V 3 , and V 5 are 
constants which depend on body shape. The Frossling 
number distribution can then be expressed as 


= 1 . 8006 ( 4 :) - 0.4237(— ) 3 + 0.0447(-) 5 
(pV)ch D u V 

(A.9) 

For the 3:1 elliptical leading edge the result is 


Fr = F 0 + F,(^) 2 + F 4 { : i ) 4 

(A.2) 

where 


Fq — — /o \/ 2 Vi 

(A- 3) 

£ 

II 

1 

GO 

to] 

(A.4) 

„ 12V 5 ( Vi \ 

4 " y/W[ V* VlVsJ 

(A.5) 

The coefficients /o,/ 2 , 04 ,h 4 depend only 
number. For Prandtl number equal to 
are [ 10 ] 

on Prandtl 
0.70 these 

fo = -0.4959 / 2 = -0.1119 

g A = -0.0977 h 4 = 0.0318 

(A. 6 ) 

For air at room temperature, it is more accurate to as- 
sume Pr = 0.71 [14, p. 269]. Rather than recalculate 
these coefficients for Pr = 0.71, it is assumed that 

„ „ / 0.71\° 4 

Fr ° 71 = Fr ° " 0 U™ J 

(A. 7) 

at each location on the body [16], resulting 
imately a one-half percent adjustment to tr 
number. 

; in approx- 
ie Frossling 


Thus, by specifying the coefficients Vj, V3, and V5 
in equation A.l, one can predict the Frossling num- 
ber distribution in the leading edge region. These co- 
efficients are, in general, dependent on overall body 
shape. Effect of overall body shape must be taken 
into consideration for an accurate prediction of the 
leading edge heat transfer. As will be shown below, 
the two body shapes considered in the present study, 
with exactly the same leading edge radius of curvat ure, 
can produce stagnation point Frossling numbers which 
differ by roughly 20%. It will also be shown that this 
difference can be completely accounted for in terms of 
the coefficient Vi for each body. From these results, it 
will also be seen that one should not expect a leading 
edge Frossling number near unity for all body shapes. 
One must know the gradient of the edge velocity at the 
leading edge, V\ to make an accurate approximation. 


- 1 . 2027 ( 4 :) - 0 . 6072 ( 4:) 3 + 0 . 2239 ( 4;) 5 

{pV)ch D D V 

(A.10) 

The edge momentum and the corresponding fifth or- 
der polynomial are shown in figure 15. The fact that 
equations A.9 and A. 10 pass directly through the nu- 
merical results indicates that a fifth order polyno- 
mial adequately represents the solution in the range 
of 0 < s/D < 1. Notice that for the 3:1 elliptical lead- 
ing edge, the gradient at the leading edge, Vi, is 33% 
lower than for the circular leading edge. 



Figure 15: Edge momentum vs streamwise distance. 

Putting these values into equations A.2-A.5, and 
making the correction for Prandtl number equal to 
0.71, equation A. 7 gives 

Fr = 0.9411 - 0.1999(-^) 2 + 0.0076(-^) 4 (A. 11) 

for the circular leading edge, and 

Fr = 0.7691 - 0.3505(|;) 2 + 0.0938(-^) 4 (A.12) 

for the 3:1 elliptical leading edge. Thus, the Frossling 
analysis predicts that the Frossling number at the stag- 
nation point will be 0.9411 for the circular leading 
edge, and 0 7691 for the 3:1 elliptical leading edge. 
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